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NONLINEAR DYNAMIC ANALYSIS OF DEPLOYING FLEXIBLE SPACE BOOMS 


by 

Paul E. McGowan and Jerrold M. Housner 


SUMMARY 


A fundamental Investigation of the planar deployment and lock-up of 
two flexible boom-type appendages which have attached tip masses and are 
connected to a central rigid body through a rotational spring is presented. 
Nondimensional parameters are identified and it is shown that, in general, 
the solution depends only on two mass ratios and one nondimensional stiff- 
ness parameter. Results are presented for boom tip deflections, deployment 
time and root moments at lock-up. A threshold value of the nondimensional 
stiffness parameter is identified beyond which boom deflections become 
large. Also, a thorough examination of the effect of nonlinear terms in the 
equations governing the deployment phase is performed. Nonlinear terms in 
the deployment equations due to kinematics and structural deformation are 
required to predict more accurately boom deflections, but retention of an 
inconsistent set of nonlinear terms leads to erroneous results. In par- 
ticular, retaining nonlinear kinematic terms while neglecting nonlinear 
structural terms can produce inaccurate results even below the threshold 
stiffness value. 


INTRODUCTION 


Many space structure designs require the on-orbit deployment of the 
entire structure or the assembly of various deployable components. The U.S. 
Space Station for example may be predominantly composed of interconnected 
deployable subassemblies. Other examples include antenna feed masts, solar 
arrays, radiators, robotic manipulator arms, and space platforms. 
Constraints on packaged volume and weight lead to very complex packaging and 
deployment schemes which require reliability in all phases of the 
deployment. Since ground testing of many of these concepts may not be prac- 
tical, dynamic analyses will be relied upon much more than in past, projects 
for verification of deployment concepts. 

Two methods proposed for deploying large space structures are exten- 
sion (or telescoping) and unfolding. Some literature exists on the analysis 
of extensions! deployment (e.g., references 1 ~ 3). However, as noted in a 

recent review of dynamic analyses for deploying space structures (reference 
4) much less attention has been given to unfolding. This is a shortcoming 
in deployment analyses because, as proposed space structures grow larger, it 
is likely that deployment will take place via unfolding or a combination of 
extension and unfolding. Thus, it is imperative that efficient, verified 
deployment analyses be developed. 





Until recently, much of the research conducted in unfolding deployment 
has been limited to cases in which the individual elements were assumed to 
be rigid during deployment (references 5-8). Reference 9, however, 
presented results of a convected finite element analysis treating the two- 
dimensional unfolding of slender, flexible beams. This analysis allowed for 
many interconnected beams and accounted for large rotations and nonlinear 
deformations. Considerable insight into the behavior of these types of 
structures was gained; however, the very general nature of the analysis 
method used makes extensive parameter studies on a computer impractical. 
Therefore, it is desirable to perform fundamental deployment studies to 
verify general analysis computer programs, to determine the effect of cou- 
pling between nonlinear kinematics and flexiblity, and to gain insight into 
the unfolding deployment process. 

The purpose of this paper is to present results of a fundamental study 
on the dynamics of deployment, examining the effects of various nonlinear 
terms arising from different analysis assumptions. The problem used to il- 
lustrate these effects is the planar symmetric deployment and subsequent 
lock-up of two flexible boom-type appendages which are elastically connected 
to a central rigid body through a rotational spring. The equations of mo- 
tion are formulated for the deployment and post-lock-up phases and general 
nondimensional parameters are identified. The deployment phase is analyzed 
with nonlinear equations, while the post-lock-up phase is analyzed with 
linear equations. Limiting cases are examined in closed form and the 
general equations are solved numerically. 

Solution accuracy as it pertains to the retention of nonlinear terms 
is also addressed. Typical results are presented for boom tip deflections, 
deployment time and peak root bending moments and comparisons are made with 
results obtained from a finite element deployment code. 


SYMBOLS 


c damper coefficient 

d vector of problem degrees-of-freedom; d = (0, v/L, w/L) 

d reduced vector of problem degrees-of- freedom when m R goes to 

m infinity; a m = (0, w/L) 

d reduced vector of problem degrees-of-freedom when boom is assumed 

1 to be rigid; d = (e, v/L) 

D damping matrix 

E Young's modulus 

f. jth approximation to the vector of nonlinear terms given by 

J eq. (4a) 

F vector of external force 
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ith component of the vector of nonlinear kinematic and structural 
forces from eq. (3;; i=1 , 2, 3, 4 

dimensional inertial coefficients which are functions of and 
m T as provided in the Appendix; 1=1,..., 6 

boom cross-sectional area moment of inertia 

nondimensional inertial terms which are I . I. , functions of m p 
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and m as provided m the Appenaix 


nondimensional stiffness parameter, k = 

dimensional rotational spring constant 

nondimensional linear stiffness matrix 

boom length 

mass of boom 

one half the root mass 

tip mass 
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mass matrix 

bending moment at boom root 

average axial load along the boom due to centrifugal acceleration 
nondimensional coordinate along the boom 

applied moment at boom root due to rotational spring and viscous 
damper 

time 

axial displacement of boom tip as shown in figure 2 
axial displacement along the boom as shown in figure 2 
vertical motion of the boom root as shown in figure 2 
transverse deflection of the boom tip as shown in figure 2 
transverse deflection along the boom as shown in figure 2 
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orthogonal axis system shown in figure 2 
rigid boom pendulum frequency. 
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nondimensional frequencies given by eqs. (12) and (13) 

uncoupled vibration frequencies given 
in eq. (6) 

U) | / 10 2 

rotational degree-of-freedom as shown in figure 2 
nondimensional time = toip 

nondimensional deployment time at which 0 =it/ 2 
value of t d for a rigid boom 


displacements parallel to the x and y axes respectively 

critical damping ratio in pendulum 
c w 


mode, x, = — 

2 / k 


ANALYSIS 


The planar unfolding deployment of two boom-type appendages which have 
attached payloads and are connected to a central rigid body through a rota- 
tional spring is depicted in fig. 1. Assuming symmetric deployment, the 
mathematical model of fig. 1b is applicable. The deployment mechanism con- 
sists of a linear prestressed rotational spring, which drives the depioyment 
and is at its neutral positon at lock-up, and a linear rotational viscous 
damper which can be used to control the motion during deployment. Thus, the 
applied moment at the root during deployment is 

T = k(n/2-0) - c0 

where k is the rotational spring constant and c is the rotational viscous 
damping constant. When the boom root has deployed through a 90 degree angle 
the root joint is assumed to lock-up and become perfectly rigid (herein 
referred to as lock-up), and the boom vibrates in the post-lock-up phase as 
depicted in fig. 1b. 

The analysis presented herein is divided into two parts, a deployment 
phase and a post-lock-up phase. The solution of the deployment phase 
provides the intial conditions for the post-lock-up phase. 


it 



Deployment Phase 


The motion and deformation of the boom , (see fig. 2), may be 
described by 0, the rotation of the boom root; v/L, the nondimensional 
translation of the boom root along the y axis; U/L, the axial displacement 
of the boom along its length; and w/L, the transverse deflection of the boom 
along its length. Also shown in fig. 2 are u, the axial displacement of the 
boom tip and w, the transverse deflection of the boom tip. 

The equations of motion describing the deployment phase (i.e., from 
rest, 6 = 0, to lock-up, e =tt / 2 ) are highly nonlinear due to the presence of 
kinematic centrifugal and Coriolis forces , the large angle rotation terms 
inherent in the mass matrix and the moderate deformations of the boom. An 
approximate solution to these equations may be found by assuming that the 
transverse boom deflection, W, and the axial boom deflection, U, have cubic 
and linear variations respectively along the boom; that is, 

2 

W(s)=ws (3~s)/2 ; U(s)=us (.1) 

where s is a nondimensional coordinate along the boom, going from zero at 
the root to one at the tip, w is the value of W ( 1 ) and u is the value of 
U(1) . 


In addition, axial deflections of the boom may be eliminated by impos- 
ing the simplifying assumption that the boom is inextensional . To first 
order nonlinear terms, this implies that the average strain over the boom 
length is zero and yields, 

2 

u/L = “(3/5) (w/L) (2) 


As a consequence of eqs. (1) and (2), the vector d containing the 
reduced system degrees-of-freedom has three entries, namely, 

1 

d = < v/L • 

(_ W/L J 

Then the equations of motion governing the deployment phase may be ex- 
pressed in the form, 

4 

Md + Dd + Kd + \ g. = F (3) 

i=1 1 

where the initial conditions on d and d are zero, and where M, D and K are 
the mass, dissipation and stiffness matrices .respectively, the sum of the 
g2 s is the vector of nonlinear kinematic and structural forces and F is the 
vector of external forces. The derivation of M, D, K, g. and F is provided 
in the Appendix. 



The matrices M, D, and K are given by, 
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are nondimensional inertial properties, 
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the percent of critical damping in the pendulum mode of the boom 
presence of the linear viscous damper at the boom root 


the pendulum frequency of a rigid boom with a rotational spring 
its root, namely, 
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The complete nonlinear force vector in eq. (3) contains many nonlinear 
kinematic and structural terms. In order to assess the importance of these 
terms the following approximations to this vectbr are introduced, 

j 

*Y - I S, ; j-1 , 2, 3, 4 (^) 

J i=1 

where each of the f. represent varying levels of approximations to the com- 
plete nonlinear fot^ce vector, and each of the g^ are given in the appendix. 
These approximation levels correspond to variousassumptions that have been 
used in the past by different investigators. 

When considering only the nonlinear terms arising from purely 
kinematic considerations of a flexible boom, the first approximate force 
vector f i may be written as, 
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(4b) 


where in this equation all terms involving products of the boom transverse 
deflection have been neglected. This order of approximation is equivalent 
to that used in ref. 10. 


When the deflection terms which were neglected in eq. (4b) are 
included, the force vector f^ takes the form 


f 1 + 


[i {(w/L) 0 +2(w/L) (w/L) ©} ] 

ww 


o 
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(4c) 


This order of approximation is the one most often used when nonlinear ef- 
fects are included; i .e . , a linear structural representation is used, often 
in the form of a summation of modes, and the resulting nonlinear kinematic 
terms are retained. 


The next level of approximation to the force vector is obtained by in- 
cluding differential stiffness along with the nonlinear kinematic terms of 
eq. (4c), thus, f^ is written as, 



( 4d) 


where the differential stiffness effect is the term in the third row of the 
vector added to f 2 . 

As a final level of approximation, when along with the nonlinear 
kinematic and differential stiffness terms, a complete set of first-order 
nonlinear structural terms are retained, the nonlinear force vector, f^, is 
defined. In accordance with eq. (4a), f^ is assembled using eqs. (A7) 
through (A10) in the Appendix. 


The effects of the various levels of approximations to the nonlinear 
force vector on the solution accuracy are examined in a later section. 
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Also from eq. (3), the vector F in dimensionless form may be expressed 
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Furthermore a dot (.) as used in eq. (3) indicates differentiation with 
respect to a dimensionless time parameter, t, which is related to time t 
as , 

i = a) t (5) 

P 

An examination of eq. (3) reveals that, for given values of root and 
tip mass to boom mass ratios, only one parameter, namely k, is required to 
determine the nondimensional vector d. 

Equation (3) may be solved by direct numerical integration using any 
of a large selection of explicit integrators. In addition, certain limiting 
cases of eq. (3) may be solved in closed form when the nonlinear force vec- 
tor is neglected, such as when eq. (3) is an uncoupled set, when the boom is 
rigid and when the root mass is very large. Each of these limiting cases is 
discussed in the following sections. 


Uncouple d Eq uations - Equation (3) represents a coupled set of 
equations . Coupl fivg'Ts __ pr ovi ded through the mass matrix and through the 
nonlinear force vector. If it is assumed that eq. (3) can be approximated 
as an uncoupled set of three equations, the nondimensional natural vibration 
frequencies with nonlinear terms neglected are, 


w =0 
v 

to = / 17k (6) 

w 

where w is associated with the pendulum mode of a rigid boom, to w is as- 
sociated with the fundamental cantilever boom frequency since the deflection 
W(s) has been chosen such that its gradient is zero at the root, and w v is 
associated with the rigid body freedom of the system in the v direction. 
The dimensional frequencies are found by multiplying the nondimensional fre- 
quencies by w . 


Rigid Boom - Equation (3) may be readily reduced to the limiting case 
of a rigid boom. This is accomplished by setting w to zero and striking the 
third row and column of M, D, and K, and the third row of F. Note that the 
nonlinear force vector vanishes when the boom is treated as rigid. The 
resulting equation is 
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As revealed by eq. (7), all nondimensional results are independent of the 
rotational spring oonstant, k. For example, as the root mass goes to 
infinity, (or equivalently v goes to zero), 6 limits to 0 , where 


0 =(ir/2) [l - 
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t J T-i 


and the nondimensional deployment time 
0=tt/2) limits to -t^, where 
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(the time to achieve lock-up, 


( 8 ) 


Infi nite Root Mass - Equation (3) may also be reduced to the limiting 
case of - an infinite root mass or equivalently to the case of a constrained 
boom root in the v direction. This is accomplished by striking the second 
row and column of M, D, and K and the second row of F. The resulting equa- 
tion is, 
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and where the damper and nonlinear terms of eq. (3) have been neglected. 

The solution of eq. (9) is 
0 = ( tt/2 ) { 1 - bcos or 1 t - (l-b)cos cd 2 t } 
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ft = 0) 1 /w 2 

Shown in fig. 3 are the variations with k of uh and w 2 . Equation (9) 
indicates that these frequencies reflect the coupling of the 0 and w de~ 
grees-of ~f reedora through the mass matrix. The effects of this coupling on 
the solution can be seen from fig. -4, which depicts the root rotation 0 ver- 
sus the nonaimens ional time, x for flexible booms with two values of the 
stiffness parameter k. Also shown in fig. 4 is the 0 history for a rigid 
boom (recall that this is independent of k). Both curves of fig. 4 for 
k=0.10 and k=0.90 exhibit the oscillatory motion of 0, however, the number 
of oscillations for the two cases are very different since the w 2 frequency 
from fig. 3 associated with the k=0.90 case is much lower than that of the 
k=0.10 case. Furthermore, as k increases, the degree of coupling between 
and w 2 increases such that the solution for 0, and hence the deployment 
time, changes. 

Provided the amplitude of the oscillations associated with w 2 are 
small, the nondimensional deployment time is essentially tt/ 2. However, as 
boom deformation increases, the amplitude of these oscillations increases 
and as can be envisioned from fig. 3> the deployment time can change 
abruptly as the boom locks up due to these oscillations. 

Referring again to fig. 3> the coupled nondimensional frequency mi 
limits to Wq as the stiffness parameter k goes to zero, (i.e., the rota- 
• tional spring stiffness vanishes or the boom becomes rigid), and limits to 

the cantilever boom frequency, / 3/k as k approaches infinity, (i.e., the 
boom root becomes cantilevered or the boom stiffness vanishes). Since aij. is 
very close to w for k less than about one, it is natural to refer to oil as 
the coupled penaulurn frequency. 


For small values of k, the coupled frequency m 2 limits to 
and for large values of k limits to w where, 


m* = [I /(I -I* ,)] 
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If the tip mass to boom mass ratio were allowed to approach infinity, to 2 
would approach the nondimensional frequency of a beam with a rotational 
spring at the root and a simple support at the tip, and if in addition, k 
were allowed to go to zero, m 2 would approach the nondimensional frequency 
of a simply supported beam. It is thus reasonable to associate w 2 with boom 
flexure. 
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Post-Lock-Up Phase 


The solution of eq.(3) at t= t d provides the initial conditions for 
the post-lock-up analysis. Whereas the deployment phase is analyzed with 
nonlinear equations, the post-deployment phase is analyzed with linear equa- 
tions with nonlinear initial conditions on W and W. It is solved herein 
using modal superposition where the modes of a beam with a guided root and 
masses at each end are used. In terms of the nondimensional parameters used 
herein, the beam equation for the post-lock-up phase is, 

w «,« + -Ck/ <1/3 + m ) ] ( W ) =0 (14) 

I ooSo i 

with boundary conditions, 

W,£g l = (1/2)m R kW(0)/[(1/3)+m T J 
W ^0) =0 


= ~ ( 1/2 )m T k W(1)/[(1/3)+m T J 

where, and m R are the ratios of the tip mass and root mass, respectively, 
to the boom mass. 


DISCUSSION OF RESULTS 


Accuracy 

The accuracy of the equations of motion (3) is considered in this 
section. Accuracy is especially important in the deployment phase which is 
analyzed with nonlinear equations. The nonlinear terms arise from boom 
flexure and enter the equations of motion through both kinematic and struc- 
tural considerations. It is often difficult for the analyst to decide which 
nonlinear terms should be retained in the analysis so as to generate a con- 
sistent formulation which produces accurate results. When, as is often 
desirable, the nonlinear equations of motion are linearized about a given 
state in order to assess stability, an inconsistent set of nonlinear terms 
retained in the equations can lead to erroneous results. 

Accuracy of deployment time, ■[_, (which is a measure of the accuracy 
in computing 6), and transverse deflection of the boom tip at lock-up, w/L, 
are examined for the case of the root mass to boom mass parameter, m R , ap- 
proaching infinity and the tip mass to boom mass parameter, m^. , equal to 
1.0. It is anticipated that the general conclusions for this case are also 
valid for other values of m u and m_. 

Discussion of accuracy is divided into two parts. First, in figs. 5a 
and 6a the solutions given by the linearized equations of motion (nonlinear 
force vector neglected) and the equations of motion with the nonlinear force 
vector (given by f^) employed are compared to the exact solution taken as 
the converged finite element results of ref. 9. Second, various levels of 



approximations to the nonlinear solution as given by through f^ are com- 
pared in figs. 5b and 6b. 


Comparison of Linear and Present Nonlinear Solutions to Exact 
Solution - The exact, linear and present nonlinear solutions are referred to 
and discussed as follows: 

(1) Exact nonlinear kinematics/ exact nonlinear structure 

(2) Linear kinematics/linear structure 

(3) Exact nonlinear kinematics/approximate nonlinear structure, 

Exact nonlinear kinematics/exact nonlinear structure - The converged results 
from the general finite element nonlinear deployment code of ref. 9 are 
taken to be exact. This program allows for unlimited deformations and rota- 
tions with the assumptions herein relaxed. 

Figure 5a displays the variation in the nondimensional deployment 
time, Tp , with the nondimensional stiffness parameter, k. Increasing values 
of k can be interpreted as either increasing the value of the rotational 
spring constant or as decreasing the bending stiffness of the boom. 
Similarly, fig. 6a displays the variation in the boom tip deflection at 
lock-up, w/L, with k. 

For small values of k the deployment time is approximately that of the 
rigid boom case, u =^/2; however, as k increases a threshold value is 
reached, near one, at which the deployment time decreases rapidly. As k ap- 
proaches infinity the deployment time goes to zero. In addition, as k 
increases beyond the threshold value the boom tip deflection of fig. 6a also 
increases and limits to L as k approaches inf inity; i .e . , no axial extension 
of the boom occurs. Sketches of the boom's deformed shape at lock-up are 
provided in fig. '6a to allow the reader to visualize what occurs to the 
physical boom as k is increased. 


L inear kinematics/linear structure - The solution for 0 and w when only 
linear terms are retained is given by eqs. (10) and (11). The deployment 
time t is found from the transcendental equation which results from setting 
6 to ir/2 in eq. (10). Also shown in figs. 5a and 6a is the variation of 
and w/L with the nondimensional stiffness for this solution. 

The value of depends on tuj , associated with pendulum motion of the 
boom, and m 2 , associated with flexural vibration of the boom. Recall from 
fig. 3, that at small values of k, w 2 is much higher than idj and thus does 
not couple with to affect deployment time. Consequently, the linear 
solution from fig. 5a is very close to the rigid boom pendulum deployment 
time (i when k is small, as was the exact solution. Furthermore, as 

shown in rig. 6a, the boom tip deflection magnitude at lock-up is small 
provided k is small. As the value of k increases, however, the two modes 
associated with au and w 2 couple together. This coupling appears to produce 
the threshold dropoff noted with the exact solution. For k above the 
threshold value of about one, the deployment time decreases rapidly to a 
value of about 0.2, (the linear equations do not permit the curve to go to 
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zero aa k goes to infinity as in the case of the exact solution), and the 
tip deflection magnitude as seen in figure 6a gets large; limiting to about 
1.6L as k goes to infinity compared to L for the exact solution. Thus the 
linear equations are in error for predicting tip deflection for very large 
values of k. 


Exact nonlinear kinematics/ approximate nonlinear structure, f ) | - When, along 
with linear terms, the nonlinear inertial, differential stiffness and first- 
order nonlinear structural terms are retained, the nonlinear vector as given 
by f^ is employed. These nonlinear terms are derived in the Appendix by al- 
lowing axial motion of the boom in such a way as to constrain the average 
axial strain in the boom to vanish to first-order nonlinear terms. 

The predicted deployment time and tip deflection magnitude for this 
solution are also given in figs. 5a and 6a respectively. The results indi- 
cate that for very small k the inclusion of the nonlinear terms has little 
effect and results are close to those of the rigid boom and linear equations 
solution. For deployment time, this solution also indicates a threshold k 
value of 1.0 and follows the linear and exact solutions closely. For tip 
deflection at lock-up, the nonlinear solution approaches a deflection of 
about 1.25L as k becomes large. This implies that there is still some inac- 
curacy in computing the boom tip deflection at very large k even though the 
boom axial strain has been constrained to first-order nonlinear terms, 
however, the error is considerably less than that given by the linear 
solution. 

In light of these results it appears that a threshold value of k ex- 
ists beyond which considerable deformation of the boom results and that the 
linear solution predicts the existence of this threshold. For the case of 
m^ very large and m equal to one, the threshold value is about one. It is 
expected that results for other mass ratios will exhibit the same character 
with different threshold values. 


Various Approximations to Nonlinear Force Vector, f ^ - Consideration 
is now given to how the various levels of approximations to the nonlinear 
force vector, eqs. (4a)-(4d), affect the solution accuracy. For comparison, 
the results from figs. 5a and 6a for the complete nonlinear force vector 
(exact nonlinear kinematics/approximate nonlinear structure, fj.) are 
repeated on figs. 5b and 6b. Each of the approximations are ref erred to and 
discussed as follows: 

(1) Approximate nonlinear kinematics/ linear structure, f 1 

(2) Exact nonlinear kinematics/ linear structure, f 

(3) Exact nonlinear kinematics/linear structure with differential 
stiffness, f 

Approximate nonlinear kinematics/linear structure,^ - When only the non- 
linear inertial terms arising from purely kinematic considerations are 
retained, the nonlinear force vector is given by f . . Some references (e.g. 
ref. 10) retain only terms of this type by neglecting products of the defor- 
mation (w/L), which is assumed small. Curves for deployment time and tip 
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deflection magnitude at lock-up with the f vector employed are given in 
figs. 5b and 6b respectively. 

The results indicate that this level of approximation predicts a con- 
siderably lower threshold value of k=0.02 than the linear and nonlinear 
solutions from fig. 5a. Furthermore as the threshold value of k=0.02 is ex- 
ceeded, the tip deflection prediction grows without bounds. 


Exac t nonlinear k inematics/linear structure, f 0 - When the kinematic terms 
which" were neglected in obtaining f^ are retained, the force vector is given 
by fp. Note that the deployment time prediction from fig. 5b for this case 
peaks at a large value after the threshold value of k = 1.0 is exceeded and 
then rather than dropping to zero levels off to a value higher than the 
rigid deployment time. In addition, the tip deflection prediction from fig. 
6b grows without bounds as the threshold value of k=1.0 is exceeded. Thus, 
this level of approximation is in poor agreement with the exact results. 

The unbounded growth in the tip deflection for the f^ and f 2 cases may 
be traced to the presence of the term in the third row of both vectors. If 
the boom were rotating at a constant velocity, this term would enter the 
stiffness matrix of the system as a negative linear term on the diagonal; 
thus subtracting from the boom bending stiffness. Clearly, for a suffi- 
ciently large angular rotation rate, boom deformation would grow without 
bound. 


This physically impossible result occurs due to the neglect of boom 
stiffening from centrifugal forces which become significant when the angular 
rate gets large. Examination of this stiffening term shows it to be of the 
same order of magnitude as the inertial term which is usually retained by 
analysts and consequently the stiffening term should also be retained. 


Exa ct nonlinear kinematics/linear structure with differential stiffness, f ^ 
When,' in addition to the exact nonlinear kinematic terms of f 2 , differential 
stiffness is included, the nonlinear force vector is given by f^. From eq. 
(^d) the coefficient of the term in the third row may be shown from the 
definitions in the Appendix to always be positive. Thus inclusion of dif- 
ferential stiffness prevents the deformations from growing unboundedly as 
shown in fig. 6b. The tip deflection magnitude is seen to limit to about 
1 . 2L as k becomes very large. Although there is still some inaccuracy in 
computing the boom deflection at the upper limit of k, the error has been 
significantly reduced by the addition of the differential stiffness term. 

Also, for k above the threshold value of 1.0 the deployment time from 
fig. 5b increases somewhat above the rigid boom solution rather than drop- 
ping off rapidly as did the linear solution. Therefore, the retention of 
the differential stiffness term in the w equation improves the prediction of 
boom tip deflection, yet it still does not yield an accurate prediction of 
deployment time, which is associated with the 0 equation. 

From the foregoing, it thus appears that retaining only some of the 
nonlinear terms arising in the equations of motion can lead to inconsistent 
results. This is especially true for the two cases (f-j and f when only 



the nonlinear inertial terms arising strictly from kinematic considerations 
are retained. 


Deployment Time 

As displayed in fig. 5a, the nondimensional deployment time of a 
flexible boom, whose value of k is less than the threshold value, is closely 
approximated by that of a rigid boom. This is generally true independent of 
whether or not nonlinear terms are retained. Although threshold values for 
various mass ratios are not provided herein, it is expected that in many ap- 
plications typical values of k will be quite small and less than the 
appropriate threshold value. For example, if a rotational spring constant 
were sized to rotate the upper and lower arm booms of the Space Shuttle RMS 
(ref. 11) through a 90 degree rotation in 15 seconds while carrying a 
32,000 pound payload, the value of k would be 0.10 or one order of magnitude 
below the threshold value found in figs. 5 and 6. Thus rigid boom deploy- 
ment times should provide reasonably accurate approximations for typical 
flexible booms. 

Figure 7 displays the deployment times for rigid booms with various 
values of the two mass ratios. From eq. (7), the dimensionless deployment 
time depends only on the terms I and I y which in turn depend only on 
the two mass ratios. The actual depIoyment V £ime may be computed via eq. (5) 
by dividing by the pendulum frequency u) p . 

Effect Of Root Damping 

The effects of a rotational viscous damper in the deployment mechanism 
for controlling rotation rate during deployment is examined in this section. 
Displayed in fig. 8 is the root rotation history during the deployment phase 
for the case of k=0.l6 and for the same mass ratios used in fig. 5. The 
curve denoted as the rigid assumption is the solution to eq. (8) and repre- 
sents the rotation of a rigid boom with 2 % damping applied. The solid curve 
provides the undamped rotation of the boom root accounting for flexibility. 
This solution is seen to oscillate throughout the deployment due to the boom 
flexural vibrations associated with the u> 2 frequency from eq. (13)- The 
dashed curve represents the root rotation for a flexible boom with 2 % of 
critical damping added at the boom root. Note that the root damper at- 
tenuates the boom flexural vibrations. This effect is due to the strong 
coupling between rigid body rotation and boom flexure. 


Post Lock-Up Root Moments 

Of interest in the post-lock-up phase is the peak root bending moment 
following lock-up. Shown in the solid curve of fig. 9 is the variation of 
the peak root moment M^ normalized by EI/L with the nondimensional stiffness 
parameter k for the same mass ratios used in fig. 5. For these results the 
range of k was limited to be below the threshold value where the linear and 
nonlinear solutions are in good agreement and the boom deflections are 
small . 
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Also shown by the dashed curve in fig. 9, is the peak lock-up moments 
under the assumption that the boom remains perfectly rigid during the 
deployment phase, and then is allowed to become flexible in the post lock-up 
phase. The trend of the dashed curve shows that the root bending moment is 
proportional to the square root of the rotational spring constant. The dif- 
ference in the two curves of fig. 9 indicates that neglect of boom 
flexibility during deployment leads to 'nonconservative predictions of root 
lock-up moments for essentially all values of the stiffness parameter. 


CONCLUDING REMARKS 


A fundamental investigation of the planar deployment and lock-up of 
two flexible boom-type appendages which have attached tip masses and are 
connected to a central rigid body through a rotational spring has been 
presented. The equations of motion were formulated for the nonlinear 
deployment phase and the linear post-lock-up phase. Nondimensional 
parameters were identified and it was shown that the solution depends only 
on two mass ratios and one nondimensional stiffness parameter. 

Accuracy of the deployment equations was examined through comparison 
with converged results from a finite element deployment code and a threshold 
value of the nondimensional stiffness parameter was identified beyond which 
the boom undergoes large flexural deformations. The linear solution 
produces accurate deployment times, but beyond the threshold does not 
produce accurate boom deflections. 

For more accurate boom deflection predictions above the threshold 
value, nonlinear effects need to be included. However, it was shown that 
great care must be exhibited in retaining nonlinear terms, for an inconsis- 
tent choice led to erroneous results even below the threshold value. In 
particular, the retention of only nonlinear inertial terms derived from 
purely kinematic considerations can lead to erroneous results. 
Unfortunately it is only these nonlinear terms which are retained in many 
analyses. Retention of a complete set of nonlinear terms as derived herein 
were generally in good agreement with converged results from a general 
finite element deployment program. 

Typical booms have small stiffness parameter values and it was found 
that for booms with stiffness parameters less than the appropriate threshold 
value, deployment times are accurately predicted under the assumption of a 
rigid boom deployment. Deployment times for rigid booms with various mass 
ratios were presented. 

It was found that a damper used to control rotational rate during 
deployment also suppresses boom flexural vibrations during deployment. This 
effect is due to the strong coupling between rigid body rotation and 
flexural deformations. 

Results were also presented for post-lock-up peak root bending 
moments. While the assumption of a rigid boom during the deployment phase 
led to accurate prediction of deployment times for small stiffness 
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parameters, the same assumption made in conjunction with allowing 
flexibility following lock-up led to underpredicting the peak root moments. 



APPENDIX 


The equations of motion may be derived using the principle of virtual work, 


6r B +6r s +6r I +6r D =o 


(A1 ) 


where, Tg, r„, r, and r D are the energies 

associated with tfte boom, the rotational spring, the inertial forces and the 
damping forces respectively, and <5 is the variational operator. 

From ref. 12, the contribution of the internal boom forces to the total vir- 
tual work may be expressed as, 


sr 


B 


. W + W ]SW 

L ,ss J ,ssss J 


+P6U }ds+{[ f W + (4) W _]6W} 


, SSS' 


3=1 


(A2) 


where, P represents the average axial load in the boom. 


Substituting the assumed deformation shapes of eq. (1) into eq. (A2) yields, 
5r D =C^r w + ~ w]6w + P6u 

B L 5L 

The contribution of the rotational spring to the total virtual work is, 

6r g -(| -0)k60 

and that of the damper is, 


6r B =cu>p060 

The contribution of the inertial forces to the total virtual work is, 

. i 


<s r _ 


Vp 


oj Ln6n + 


I 

where, 

5 = [sL+U(s)]sine - W(s)cos0 
n = v + [sL+U(s)]cos0 + W(s)sin0 


5S£]ds + m T [n6n + 5SO s=1 


(A3) 

(A4) 

(A5) 


After differentiating eqs. (A^) and (A5) with respect to dimensionless time 
x, utilizing the assumed shape functions of eq. (1) and substituting the 
results into eq. ( A3 ) » the virtual work of the inertial forces can be 
rewritten as, 

«r I /(m^«p)- { CH 1 L +H 2 w 2 + H 1 uL(2+u/L)]0 

+ (H c wcos0-H ll Lsin0)v + 2H.L( 1 +u/L)u0 



. • 2 " '* « • • 2 
+ { -H 1 L 6 + H^we + u + 2H 3 W0 -Hjfl u 

+ H^VCO30}6U 

+ { -H„L( 1 +u/L) 0 +H„w + H vsin0 
3 - d o 

• 2 • • l 

- H^W0 - 2H 3 ue}5w 

+ { (H^-wcose-H^sine) 0 + H^wsinO + H^v 

• 2 • • 

+ (H^cose-H^wsine ) 0 + 2H^ewcos0 

• 2 • • I 

+ H i) [u(esine+0 cose) + 2eusine + ucos 0 | 6 v 

where, 

H^I/3 + m T ; H 2 =33/1 i <0 + try H^II/HO + m T ; 

H JL( =1 /2 + m T J H 5 =3/8 + m T ; H & =1 + m T + m R 

Substituting the virtual work contributions into eq. (A1) and setting the 
resulting coefficients of the variational terms, 60, 6u, 6w and 6v, to zero, 
yields the equations of motion. The equations of motion involve the de- 
gr ees-of-f reedom 0, w, v, u and the boom tensile load P. The equation of 
motion associated with axial motion, 6u, provides the value of p as, 

_ 2 # 2 *• .. 2 

P=rn B° J p^ H 1 Lj0 ~ H 3 W0 “ 2H 3 W ® -H i u +H 1® u ~H ^VGOSQ ) (A6) 

where the right hand side of eq. (A6) is recognized as an approximation to 
the centrifugal force in the boom. 

A simplifying assumption is now made that the boom is inextensional . 
To first order nonlinear terms this implies the relationship of eq. (2). 

As a consequence of eqs. (A6) and (2), the equations of motion a3 
given by eq. (3) are produced where the terms of the system mass matrix are 
given by, 


W 8 , ; I vv‘ H 6 /H , 1 ! 

W 1 *, ! I ev- H i. /K i 

Grouping the nonlinear terms in accordance with eq. (4a), the g i vec- 
tors can be written as 





V J 

where the nonzero components are 
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g.^= I (w/L) (v/L)cose 
1 vw 

(v) ** • 2 •• 

g„ =1 [6(w/L)cos0 -0 (w/L)sin0 + 20(w/L)cos0] 

1 vw 

-I 0y 0 GOS0 

( w ) • 2 

So G) “ I ,J(w/L) 2 0 +2(w/L) (w/L) 0] 

2 WW 

(w) . 2 

g^” = ( 6/5) (w/L) 0 


(A7) 

(A8) 

(A9) 


gj 0) = -(6/5) (w/L) 2 0[1-(6/5 )(w/L) 2 ] 

-(3/5)I 0w (w/L)[(w/L) (w/L) +2(w/L) ] 

-(12/5) (w/L) (w/L) 6[1- (3/5) (w/L) ] 

+(3/5)I 0y (w/L) (v/L)sin0 

(v) 2 - . 2 

= 1 0V .[ ( 3/5 ) (w/L) ( 0sin0 + 0 cos0) 

+( 1 2/5) 0(w/L) (w/L)sin0 

- ( 6/5 ) L { (w/L) (w/L) +(w/L) cos0}] 

gj = (18/25) (w/L) [2(w/L)-0 (w/L)+2(w/L) /(w/L)] 

-(6/5)1. (v/L) (w/L)cos0 

0V 

-(3/5)I Qw 0(w/L) 2 (A10) 

where the superscripts (0), (v) and (w) refer to the degrees-of-freedom 
which the components of the g^ vectors are associated with. 
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Fig. 1 Fundamental deployment problem. 








Fig. 3 Variation of nondimens ional frequencies 
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Fig. 4 Root rotation history during deployment. 
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Fig. 5 Variation of nondimensional deployment time with nondimensional stiffness. 
Fig. 5(a) Comparison of selected approximations with exact solution. 
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Fig. 5(b) Effect on predicted solutions to neglecting various orders of 
nonlinear kinematics and structural terms. 
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Fig. 6 Variation of boom tip deflection at lock-up with nondimensional stiffness. 
Fig. 6(a) Comparison of selected approximation with exact solution. 
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Fig. 6(b) Effect on predicted solutions to neglecting various orders of 
nonlinear kinematic and structural terms. 
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Fig. 7 Nondimensional deployment time versus mass ratios for rigid booms. 
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Fig. 8 Effect of root damping on boom vibrations during deployment 
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